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Abstract 

We introduce an algorithm aimed to reduce the dimensions of 
q-( Hilbert space. It is used here in order to study the behaviour of low 

energy states of strongly interacting quantum many-body systems at 
first order transitions and avoided crossings. The method is tested on 
different frustrated quantum spin ladders with two legs. The role and 
importance of symmetries are investigated by using different bases of 
states. 
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1 Introduction. 

The investigation of the spectral properties of strongly interacting micro- 
scopic many-body quantum systems necessitates the use of non-perturbative 
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approaches. Most of these rest on the renormalization group concept in- 
troduced by Wilson [T] and universally used since in all fields of quantum 
physics. This is the case for the study of lattice quantum systems for which 
the Real Space Renormalization Group (RSRG) El El and the Density 
Matrix Renormalization Group (DMRG) j3J Ej have been introduced in 
order to reduce the dimensions of the systems. More generally the particular 
property of the renormalization methods is the fact that they can reveal the 
existence and location of phase transitions which are characterized as fixed 
points, i.e. the coupling strengths which define the Hamiltonian or the ac- 
tion of the systems and evolve during the renormalization process (running 
coupling constants) stay constant at these specific points, both in the case 
of discontinuous and continuous transitions |S] . 

The reduction process has mainly been applied in d- dimensional real and 
momentum space. In the present approach we introduce a reduction algo- 
rithm which operates in Od Hilbert space. The spectral properties of quantum 
systems are obtained through the diagonalization of a many-body Hamilto- 
nian acting in a complete, in general infinite or at least very large set of basis 
states although the information of interest is restricted to the knowledge of 
a few low-energy states. 

In order to avoid the diagonalization of very large matrices we recently 
proposed a non-perturbative approach which tackles this question 0. We 
implemented it in the study of strongly interacting systems like frustrated 
two- leg ladders [TU]. The approach consists of a step by step reduction of 
the size of Hilbert space by means of a projection technique. It induces 
a renormalization process in the spirit of former work ^T] ^1 ^]. Its 
advantage over other reduction procedures lies in the fact that it is applicable 
to all types of microscopic quantum systems since it works in Hilbert space 
like in the procedure developed in ref. [T3] . 

In the present work we concentrate on the use of the algorithm in order 
to characterise and study interacting systems in the vicinity of discontinuous 
(first order) phase transitions and avoided crossings which are signalized by 
the existence of fixed points at which coupling constants stay constant during 
the reduction process. We investigate the ability of the algorithm to perform 
the reduction of the Hilbert space dimensions at the location of these fixed 
points and to signalize the existence of phase transition points. 

The outline of the paper is the following. In section 2 we recall briefly 
the main lines of the space reduction algorithm which has been developed 
elsewhere jH] and show the implication of the existence of fixed points on the 
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coupling strengths which enter the Hamiltonian of the system. In section 3 
we test the algorithm in the neighbourhood of fixed points which correspond 
to discontinuous transitions and avoided crossings in frustrated spin ladders. 
Conclusions are drawn in section 4. 



2 Formal algorithm and determination of fixed 
points 

2.1 Space reduction procedure. 
2.1.1 General concept 

We consider a system described by a Hamiltonian 

(g[ N ^ , , ■■■,g < p N ^ which depends on p coupling strengths 

^g[ N \gi N \''' > 9p N ^ l—> 9^^ an d acts in a Hilbert space HS N ^ of dimen- 
sion N. H( N > has N eigenvalues \Xi{9 )> i = 1, • • • , -N} and eigenvectors 

^(rt.^v-.iv}. 

If the relevant properties of the system are essentially located in a sub- 
space 7i( M ) of Ti.^ (M < N) it makes sense to try to define a new effective 
Hamiltonian H^ M '(g^ M ') whose eigenvalues reproduce the selected subspace 
and verifies 

HM(gM)\HfW(gM)) = Hg^)\^ M \g^)) 
with the constraints 

H9 {M) ) = \{9 {N) ) (1) 

for i = 1, M. If this can be realized Eq. (|1J) implies a relation between the 
coupling constants in the original and reduced space 

n (M) _ , f (N) IN) (jv)^ 

9k — Jk{9i ,92 ) •••># 1 
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with k = 1, ...,p. We show below how this effective Hamiltonian H^ M \g^) 
can be constructed. 
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2.1.2 Reduction algorithm and renormalization of the coupling 
strengths. 

Consider a system described by a Hamiltonian depending on a unique cou- 
pling strength g which can be written as a sum of two terms 

H = H + H 1 (g) (2) 

The Hilbert space TiS N ' of dimension N is spanned by a set of basis 
states i — 1, • • • ,N} which may be chosen as eigenstates of H . An 

eigenvector l^i ) can De decomposed on this basis 

JV 
i=l 

Using the Feshbach projection method jT^j Tv- N ' is decomposed into sub- 
spaces by means of the projection operators P and Q, 

H (N) = p H (N) + Qn (N) 

Here the subspace P7v- N ^ is chosen to be of dimension dim PTC^ N ' = N — l 
by elimination of one basis state. The projected eigenvector P|\1/[ 7V ' ) ) obeys 
the Schroedinger equation 

H tf) (^>)P\< > ) = A< J V|*f , > • (4) 

where H e ff(X^) operates in the subspace PTi.( N \ It is a nonlinear function 
of the eigenvalue A^ [2J, the eigenenergy being equal to the eigenvalue 
associated to l^f^) in the initial space TiS N \ In practice the coupling 
which characterizes the Hamiltonian in 1iS N > is aimed to be changed in 
such a way that the eigenvalue in the new space T^^ -1 ) is the same as the 
one in the complete space 

\[ N ~ 1] = \[ N) . (5) 

Hence the reduction of the vector space from iV to iV — 1 results in a 
renormalization of the coupling constant from g( N > to g^ N ~^ preserving the 
physical eigenenergy A^ , i.e. A^ = A[ 7V_1 ' ) = Ai. The determination of 
g(N-i) m eans of the constraint expressed by Eq. (j3J) is the central point 
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of the procedure. If the dependence on g is linear, i. e. H = H + gHi, 
it is obtained as a solution of an algebraic equation of the second degree in 
g^-V, see Eq.(17) in Ref. 0. 

Starting from H^ -1 ) the reduction procedure is iterated step by step in 

decreasing dimensions of the Hilbert space, N ^ N — 1 N — 2 m At 

each step g is renormalized to a new value. In the limit where the dimensions 
of 7i are very large one may go over to continuous space dimensions, (n, n — 
1) I— > (x, x — dx), and the evolution of g will be given by a flow equation, 
see Eq.(24) in Ref. jU| and Ref. ^H]. The algorithm can be generalized to 
Hamiltonians depending on several coupling constants. 

2.1.3 Remarks. 

• The implementation of the reduction procedure asks for the knowledge 
of Ai and the corresponding eigenvector ) at any size k of the vector 
space. The eigenvalue is in principle fixed as being the physical ground 
state energy of the system. We use the Lanczos algorithm which allows 
to determine Ai and |^^) I2HJ 1211 ■ Consequently this algorithm 
has been implemented at each step of the reduction process [TO] . 

• The process does not guarantee a rigorous stability of the eigenvalue 
Ai. Indeed one notices that which is the eigenvector in the 
space H^ 1 ^ and the projected state P|\I>^) of |^^) into 7-^ fc_1 )may 
differ from each other. As a consequence it may not be possible to keep 

^ rigorously equal to A^ = Ai. In practice the degree of accuracy 
depends on the relative size of the eliminated amplitudes a^(g^). We 
shall come back to this point when the algorithm will be implemented 
in numerical tests. 

• The reduction procedure needs a fixed ordering of the sequentially elim- 
inated basis states. This ordering may be chosen by following different 
criteria. Here the states are arranged according to increasing energies 
£j = (^\ N ^\H\^ N ^) and eliminated starting from the one which corre- 
sponds to the highest energy at each step of the procedure. 

• The algorithm will be used in applications of the procedure on explicit 
models, here frustrated spin ladders, and tested in different symmetry 
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schemes. As already mentioned we concentrate on aspects related to 
systems close to phase transition points. 



2.2 Fixed points 

The eigenvalues Xk(g) of H(g) = H + Hi(g) are analytic functions of g which 
may show algebraic singularities [T71 HS1 12H| at so called exceptional points 
g = g e . Exceptional points are first order branch points in the complex g 
- plane which appear when two (or more) eigenvalues get degenerate. This 
can happen if g takes values such that = t\ where = ($>k\H\$>)~). In a 
finite Hilbert space the degeneracy appears as an avoided crossing for real 
g. If an energy level belonging to the PTC subspace defined above crosses 
an energy level e/ lying in the complementary QTi subspace the perturbation 
development constructed from H e ff(E) diverges [20j . 

Due to symmetry properties physical states can get degenerate in energy 
for real values of g e . They correspond to first order phase transitions. 

Exceptional points are defined as the solutions of [THJ 120] 

f(X(g e )) = det[H(g e ) - X(g e )I] = (6) 

and 

^ U=A( Se ) = (7) 

It is possible to show that exceptional points are fixed points of the cou- 
pling strength g which stay constant during the space reduction process. 
If {Aj(g)} are the set of eigenvalues the secular equation can be written 

as 

N 

n(A-A,) = 0. (8) 

i=l 

Consider A = A p which satisfies Eq. (0). Eq. (|7J) can only be satisfied if 
there exists another eigenvalue \ q = \ p , hence if a degeneracy appears in the 
spectrum. This is the case at an exceptional point. 

(k) 

If the eigenvalue A^- , k — N 7 N — 1, ... which is either constant or con- 
strained to take the fixed value Xj gets degenerate with some other eigenvalue 
X\ k \g = g e ) in the space reduction process this eigenvalue must obey 

= A?V.) ( 9 ) 
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which is realised in any projected subspace of size k and I containing states 
and Going over to the continuum limit for large values of iV as 

introduced above and considering the subspaces of dimension x and x + dx 

\(g e (x),x) = (Vi(g e (x),x)\H(g e (x))\*i(g e (x),x)) 

verifies 

dX i _ Q _ d\j{g e {x),x) 
dec doc 

Consequently 

d\i d\idg e 
dx dg e dx 

Since Xi(g e (x),x) does not move with x in the space dimension interval 
(x, x + dx) the first term in Eq. (|TT|) is equal to zero. Hence in general 

^#0 and ^ = (12) 
og e dx 

The present result is general, it is valid for both continuous and discon- 
tinuous quantum phase transitions. Avoided crossings reflect the existence 
of continuous transitions in the limit of systems of infinite size. In the case 
of discontinuous (first order) transitions the energies of the physical states 
(l^i)} show a real degeneracy point, i.e. their energies cross at real values 
of g = g e as already mentioned above. 

In the following we shall consider Hamiltonians linear in g, H(g) = H + 
gHi. A sufficient condition for possible level crossings is given by 



[#o,#i] = (13) 

i.e. H and Hi can be simultaneously diagonalized jHHH]. In this case if H 
is diagonal in the {\&i), i = 1, • • • , N} basis of states it is also an eigenbasis 
for H and 



d\j(g e (x),x) _ d(^j\H + g^H^j) 
dx dx 
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Since H and Hi do not depend on x, dXi/dx = implies dg e /dx = 
for this specific form of the Hamiltonian which is consistent with the general 
result Eq. (IT2I) . 



2.3 Remarks. 

• In the general case H(g) = H + Hi(g) g(x) = g e evolves through a 
flow equation which is more complicated than Eq. (24) in f*|. We 
shall restrict our numerical investigations to Hamiltonians which show 
a linear dependence on g. 

• The energy level crossings can occur for any two energy eigenstates of 
the spectrum. We shall consider level crossings or avoided crossings 
between the ground state and an excited state as well as the case of 
crossings or avoided crossings of excited states. 

The present formal considerations are now used in numerical applications. 
We analyze the behaviour of the reduction algorithm in the neighbourhood 
of fixed points and show to what extent it allows to detect fixed points in an 
application to two-legged frustrated spin ladders for different choices of the 
basis of states. 

3 Application to quantum transitions in frus- 
trated two-leg quantum spin ladders. 

3.1 The model 

3.1.1 SU(2)-symmetry framework. 

Consider spin-1/2 ladders |2"T) |2"2"] shown in Fig. 1 and described by Hamil- 
tonians of the following type 



L 



Si^Si 2 -\- J i y S 



i=l 



<ij> 



<ij> 




(15) 
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li 2i 3i 4i 5i ■•■ l! s /l= l/2 L 




1 2 3 4 5 / L 

Figure 1: Top: the original spin ladder. The coupling strengths are indicated 
as given in the text. Bottom: The ladder in the SO(4) representation. See 
the text. 



Working in a representation with fixed total magnetic projection M tot the 
basis of states is spanned by the vectors {|$fc), k — 1, • • • , N} where 

2L 

\<S> k ) = 1 1/2 mi,. ..,1/2 m^. ..,1/2 m 2L ^m i = M tot ) 

i=i 

and {m^ = +1/2,-1/2}. 

The indices 1 or 2 label the spin 1/2 vector operators Si k acting on the 
sites i on both ends of a rung, in the second and third term i and j label 
nearest neighbours, here j — i + 1 along the legs of the ladder. The fourth 
and fifth term correspond to diagonal interactions between nearest sites lo- 
cated on different legs. 2L is the number of sites on a ladder. Here we fix 
Jic = Jic = Jc- The coupling strengths J t , J/, J c are positive. In the sequel 
we restrict our analysis to the case where M tot = 0. 

In the present applications the renormalization is restricted to a unique 
coupling strength, see Eq. (J2J). It is implemented here by putting H = 
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and HW = £ (A °#i where gW = J t and 



L 

i=l <ij> 

+lc^2(s h s j2 + s h s h ) (16) 



<ij> 



where 7« = Jj/Jt, 7 C = J c /Jt- These quantities are kept fixed and g^ N > = J t 
will be subject to renormalization in the reduction process. 



3.1.2 SO(4)-symmetry framework. 

The basis of states may be written in an SO (4)-symmetry scheme. By means 
of a spin rotation [221 

Si^^Si + Ri) . (17) 

s l2 = \{S l -Ri). (18) 
the Hamiltonian Eq.(15) can be expressed in the form 



H ^ = l^-^+^E^- ( 19 ) 

i=l <ij> 



+J2 RiRj 



<ij> 

which reduces the ladder formally to a chain, see bottom of Fig. 1. Here 
J\ = (Ji + Jc)/2, ^2 = (^ — ^ c )/2, Jic = ^2 C = <^c- The components 
+S) ,s!> , and r!- ,r\ , R^ of the vector operators Si and Ri are the 
50(4) group generators and < ij > denotes nearest neighbour rung indices 
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q (z) _ y (ll)(ll) v (l-l)(l-l) 



itf } = V2(X^ 1)m - ^ 00)(1 - 1) ) = ^ 
R {z) = _(x (10){00) + x (00){10) ) 

where 

x^ s w = \sm)(s'M\ 

with states defined as 

\SiMi) = ^ (V2 mi 1/2 m^M^l^ m^ll/2 m 2 ), 

mi,m2 

along a rung are coupled to Si — or Si — 1. Spectra are constructed in 
this representation as well as in the 577(2) representation. 

In a representation with fixed total magnetic projection M tot = the 
basis of states is then spanned by the vectors 

L 

i=l 

and {Mi = 0, +1, — 1} 

In the sequel we analyse the behaviour of the system in both symmetry 
frames. 

3.2 Fixed points in the SU(2)-symmetry framework. 

3.2.1 First order phase transitions of the two-legged spin ladder. 

At first order transitions which happen at level crossings for real g the ampli- 
tudes of the wavefunctions are expected to acquire weights of the same order 
of magnitude over a large number of basis states. The analyses performed 
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in ^U] show that the use of the algorithm at fixed points should work as a 
stringent test of the method. 

The first application concerns a system described by the Hamiltonian 
ff( s i s ) given by Eq. (fT5j) using an S'f/(2)-symmetry basis of states given be- 
low Eq. (|15p. A crossing between a rung dimer phase and a Haldane phase 
appears for Ji = J c when J t / Ji — 1.401 in the case of an asymptotically large 
system j^D I21I2HI H3 EDI EI] • The ratio depends on the size of the sys- 
tem. The existence of a first order transition is analysed below, the coupling 
constant g = J t is expected to stay constant at the level crossing point. 

3.2.2 Application of the reduction algorithm at fixed points. 

We consider ladders with L = 6. Several crossings between energy levels can 
be observed in Fig. 2(a) which shows the evolution of the energies of the four 
lowest states in the M tot = subspace as a function of g = J t . The crossing 
between the ground state energy e x and the energy of the first excited state 
ti corresponds to Jt/Ji — 1.23. 

The first test corresponds to > 6. Then the three lowest excited 
states corresponding to e 2 , e 3 , e 4 get rigorously degenerate which generates a 
continuous transition line |29j . Fixing J t = 10 Fig. 2(b) shows the evolution 
of the four lowest states as a function of the size N of the Hilbert space 
following the algorithm described in Refs. The initial space dimension 

is N = 924. The stability of the spectrum is remarkable down to N < 100. 
This stability reflects in the constancy of Jt over the same dimensional range, 
see Fig. 2(c). For lower values of N deviations appear. They are due to the 
approximations inherent to the projection method as noted in the second 
remark of subsection 2.1.3 and Refs. [5J EH]- The limit of constancy of Jt 
indicates the minimum dimension of Hilbert space in which diagonalization 
will lead to the reproduction of the low energy part of the spectrum, i. e. 
the ground and first excited states. The behaviour of the spectrum shown 
here is observed for any value of J t > 6, i. e. all along the degeneracy lines. 

Fig. 2(d) gives the evolution of the entropy s per site defined as 

1 N 

s=~j2 p * inp t with ^=i^r ) i*i N) >r=i«2°i a (2°) 

at the fixed point J t = 5 which corresponds in Fig. 2(a) to the crossing 
of the ground state e\ with the first excited state e 2 . Here {a^} are the 
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amplitudes of the components of the ground state wavefunction developed 
on the basis of states {|$j >} which span the space of dimension N. The 
step discontinuity signals the transition characterized by a strong change in 
the structure of the lowest state. The characteristic singularity observed at 
this value of J t is conserved as long as the ground state keeps stable during 
the space reduction process. 

At the exact location of the fixed point the instability of the spectrum 
is sizable and the coupling constant Jt at this crossing point stays constant 
over a smaller interval of values of N. A closer inspection shows that this 
instability might be related to a numerical difficulty in the renormalization of 
J t . Indeed the coefficients of the algebraic second order equation which fixes 
it get accidentally vanishingly small at this place and consequently lead 
to strongly unprecise values of the roots of the equation. This corresponds 
to a pathological situation which may not be significative in the general case. 
Indeed, in the close neighbourhood of the fixed point, J t stays stable over a 
much larger interval of values of N when J c = 3.8 7^ Ji (see Figs. 2 (e) and 
2(f)). 

3.3 Fixed points in the SO(4)-symmetry framework - 
First order transitions 

The reduction algorithm is next applied to the same system as above but 
described by the Hamiltonian H( S,R > given by Eq. (19) with a basis of states 
>} written in the 5*0(4) symmetry framework introduced in subsection 
3.1.2. The spectrum given in Fig. 3(a) in the M tot = subspace is the same 
as in the case of the SU(2) symmetry framework as it should be. However 
the behaviour of the numerically generated spectrum at different transition 
points is quite different. 

At the crossing point between the ground state and the first excited state 
in the M tot = subspace which occurs at J t = 5 (see Fig. 3(a)) the ground 
state energy e\ and Jt remain stable all along the space dimension reduction 
procedure as seen in Figs. 3(b) and 3(c). This remarkable stability can be 
explained by the fact that the ground state wavefunction is strongly domi- 
nated by a small number of states in the 5*0(4) basis. It is however lost for 
the first excited state with energy e 2 which moves abruptly and stays then 
again constant generating successive plateaus over more or less large intervals 
in N, Fig. 3(b). This shows that dXi/dx — 0, (i — 2, 3, 4) is indeed preserved 
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by steps, but not necessarily Aj which jumps by steps over finite intervals of 
space dimensions. The jumps in {Aj} may be related to the elimination of 
non negligible components of the wavefunction during the reduction process. 

Fig. 3(d) shows the behaviour of s the ground state entropy per site which 
behaves like in the SU (2) scheme but is quantitatively smaller. This is due 
to the fact that the wavefunction amplitudes are less equally distributed here 
than in the SU (2)-scheme as mentioned above. Figs. 3(e) and 3(f) show the 
behaviour of the spectrum and coupling constant J t in the close neighbour- 
hood of the transition point. One observes that the evolution of the energies 
is smoother than at the transition point itself and the coupling constant in- 
creases slightly with decreasing N. Some curves in the figures are drawn with 
a finite width in order to facilitate the observation of the stepwise evolution 
of the corresponding quantities. Degeneracy of the states and the consequent 
constancy for J t > 6 is observed along e 2 = ^ = &i which corresponds to 
a transition line as seen in Figs. 4(a) and 4(b). The constancy of these 
quantities is preserved over the whole range of space dimensions N, except 
for {e^j's at small N. But {e^j's stay more and more constant up to the 
smallest values of N with increasing J t , see Figs. 4(c) and 4(d). This can be 
explained by the fact that the wavefunctions gets more and more dominated 
by a small number of basis states with increasing J t . Evidently the robust- 
ness of the spectrum is stronger in the present symmetry scheme than in the 
case of SU(2). 



3.4 Application of the reduction algorithm at a con- 
tinuous transition: avoided crossings. 

As already mentioned continuous transitions reduce to avoided crossings in 
finite systems. States get degenerate at complex values of this parameter. 
Genuine transitions with real parameters cannot be explicitly seen in numer- 
ically determined spectra of finite systems. Avoided crossings are not easy to 
locate. They occur at complex exceptional points and the present numerical 
procedure is aimed to follow the evolution of real running coupling constants. 

Fig. 5(a) shows the spectrum of the ladder for specific values of the 
coupling constants. One observes several possible avoided crossings which are 
rather close to each other in the interval 4 < J t < 9. The typical behaviour 
of the spectrum is shown in Fig. 5(b) for J t = 6.6891. Fig. 5(c) gives a 
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quantitative estimate of the energy fluctuations of excited states. Here 



The spectrum and Jj are relatively stable down to iV ~ 200. Stability 
is lost below this value. The same is true at other avoided crossing points. 
These points are difficult to locate, the coupling J t is complex there and our 
procedure does not fix their imaginary part. It may be that clearer signals 
can be observed for larger systems since then the gap at crossings gets smaller 
(in principle tends to zero in case of a continuous transition) and hence leads 
to a reduced the imaginary part of the coupling constant. 

Remarks: 

• In the present calculations the system has open boundary conditions. 
For odd L the entropy shows the same characteristic discontinuity as 
observed for first order transitions in Figs. 2(d) and 3(d) at the transi- 
tion point. For L even the discontinuity goes over in a finite peak which 
recalls a continuous transition in a finite size system. This shows that 
a naive interpretation of observables in finite systems can lead to erro- 
neous interpretations of the order of a transition. 

• Expected crossings may not necessarily signal continuous transitions in 
the limit of infinite systems. In this limit the order may effectively be 
different. This phenomenon has also been observed in classical systems 
see f.i. and references therein. 

4 Summary, conclusions and outlook. 

We used an algorithm which aims to reduce the size of the Hilbert space 
of states describing strongly interacting systems. The reduction induces a 
renormalization of coupling constants which enter the Hamiltonians of the 
systems, here frustrated two-leg quantum spin ladders. The robustness of 
the algorithm has already been tested in former work [TU] . 

We applied the algorithm at the location and in the neighbourhood of 
first order transition points and lines, and in the vicinity of avoided cross- 
ings which correspond to second order transitions in infinite systems. The 
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analysis has been pursued in two different symmetry schemes. As it may be 
expected the behaviour of the spectrum and the renormalized coupling pa- 
rameter depend on the symmetry framework. Indeed, the description of the 
system depends crucially on the details of the wavefunctions of the different 
energy states and their structure is different in different representations. 

In the case of first order transitions we showed that the renormalized cou- 
pling constant J t may indeed be numerically stable over a large set of dimen- 
sions of Hilbert space and signal the presence of a transition as predicted by 
the theory considerations developed in section 2.2. Strong instabilities may 
appear due to accidental numerical pathologies as mentioned in section 3.2.2. 

The stability is larger in the case of an S'0(4)-symmetry representation 
where the number of large components is reduced. In the present parameter 
space the low-lying eigenstates are strongly dominated by a small set of 5*0(4) 
basis states. This may also explain the stronger stability of degenerate states 
at the crossing points and along transition lines. 

Avoided crossings have also been investigated. The presence of these 
points is difficult to locate. It may be due to the fact that the level cross- 
ing point occurs for a complex coupling constant which is not detected in 
the present algorithm. The precise understanding of these points related to 
numerics requires further work which lies outside the scope of the present 
investigations. 

It may be mentioned here that there exists other methods which are 
aimed to detect crossing and avoided crossing points. One of them relies on 
discontinuities in the entanglement properties of wavefunctions which 
have to be known at crossing (first order transitions) or avoided crossings 
(continuous transitions). A second approach relies on an algebraic method 
|3*B] which works very nicely in the case of small systems. It is not clear 
however whether its use can be easily applied to very large Hilbert spaces 
such as those which correspond to realistic quantum spin systems. 

In summary, the present investigations show that first order transition 
points in quantum spin systems may be correlated with strong fluctuations 
in the energies of low-lying excited states. The presence of these points 
are signalized by the constancy of the strength of the couplings which enter 
the Hamiltonian along the dimensional reduction procedure of Hilbert space. 
This is predicted by theoretical considerations and verified, at least up to the 
point at which numerical stability gets lost which may happen at different 
stages of the reduction procedure depending on the structure of the ground 
state wavefunction. 
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Figure 2: Si7(2)-symmetry scheme. The {e^, i = 1,2,3,4} are the energies 
per site of the ground and lowest excited states. N is the size of the Hilbert 
space, s the entropy of the ground state per site. The number of sites is L = 6 
along a leg, Ji = J c ~ 4.07. (b) and (c) correspond to J t = 10. Figs.(e) and 
(f) correspond to J c = 3.8 ^ J/. Broadened lines are drawn for the sake of 
readableness. See discussion in the text. 
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Figure 3: SO(4)-symmetry scheme. The {e^, i = 1,2,3,4} are the energies 
per site of the ground and lowest excited states. N is the size of the Hilbert 
space, s the entropy of the ground state per site. The number of sites is 
L = 6 along the chain, (a) - (d): Ji = J c . In (b) and (c) J t = 5 (e) - (f): 
J t — 5, Ji 7^ J c , J c = 3.8. In both cases Jt/ Ji — 1.23. Broadened lines are 
drawn for the sake of readableness. See discussion in the text. 
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Figure 4: SO(4)-symmetry scheme. The {e^, i — 1,2, 3, 4} are the energies of 
the ground and lowest excited states per site. The number of sites is L = 6 
along the chain, Ji = J c ~ 4.07. (a) and (b) correspond to J t = 6, (c) and 
(d) correspond to J t = 10. See discussion in the text. 
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Figure 5: Si7(2)-symmetry scheme. The {e^, i = 1,2,3,4} are the energies 
per site of the ground and excited states. N is the size of the Hilbert space. 
The number of sites is L = 6 sites, J t = 6.891, J[ = 5, J c = 3. See discussion 
in the text 
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